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We study dynamical mass generation and the resultant helical spin orders in topological Dirac and Weyl 
semimetals, including the edge states of quantum spin Hall insulators, the surface states of weak topological 
insulators, and the bulk materials of Weyl semimetals. In particular, the helical spin textures of Weyl semimetals 
manifest the spin-momentum locking of Weyl fermions in a visible manner. The spin-wave fluctuations of the 
helical order carry electric charge density; therefore, the spin textures can be electrically controlled in a simple 
and predictable manner. 
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Relativistic electrons governed by the Dirac equation had 
been thought to be remote from condensed matter physics. 
The developments in the last decade, especially the dis¬ 
covery of graphenelH and topological insulators (TIs)||2- 
0 ], however, have established the ubiquitousness of Dirac 
fermions in condensed matter. The topologically protected 
surface states of topological insulators are generally massless 
Dirac fermions, meanwhile, the bulk states of many topo¬ 
logical insulators can be described by massive Dirac equa¬ 
tions. More recently, massless Dirac yl-lst] (and Wevli9l 426ll ) 
fermions have also been discovered in bulk materials (a recent 
experimentally discovered material is the TaAs class ll27H35n ). 


The interactions among the nominal massless Dirac 
fermions, if sufficiently strong, can dynamically generate 
a Dirac mass and fundamentally change the properties of 
fermions. This phenomenon was first studied in the context of 
particle physics !^ . In this paper we investigate the physical 
consequences of dynamical mass generation in several classes 
of topological Dirac metals. We find that the dynamically 
generated masses manifest themselves as helical spin orders. 
These types of spin order have attracted considerable interest 
in other materials and modelAs we shall show, their 
emergence in topological Dirac semimetals is a quite robust 
phenomenon, which is independent of material details. Phys¬ 
ically, the helical spin orders result from the spin-orbit cou¬ 
pling. Their descriptions as Dirac fermions permit a unified 
treatment. 


We present the results for three classes of materials. The 
first example is the edge state of quantum spin Hall (QSH) 
insulators i43l - l45n . Dynamically generated Dirac mass man¬ 
ifests itself as helical order at the edge [see FiglT]. Quan¬ 
tum fluctuations, however, can destroy this order. Proxim¬ 
ity to other materials (e.g. proximity to certain superlattice 
structures) can stabilize the helical order. The second ex¬ 
ample is the surface states of weak topological insulators, 
which are closely related to the first example, while avoiding 
the strong quantum fluctuations because of the higher dimen¬ 
sionality. The third example is Weyl semimetal in the mag¬ 
netic field. Here, the helical order depends on the direction 
of the magnetic field in a specific manner, which provides 


a sharp experimental signature for the identification of Weyl 
semimetals lf^ . The second and the third examples will be 
our focuses. 

Helical order at QSH edge. The edge of QSH accommo¬ 
dates two counterpropagating modes, whose spin is locked 
with the propagating direction ll^ld^ . The Hamiltonian reads 

HiX) = ^ Vpik - kF)c^^^c^k - ^ Vpik + kf)c^^^Clk ( 1 ) 

k k 

where the chemical potential has been absorbed into the def¬ 
inition of kf. This form of Hik) is dictated by time reversal 
symmetry. We can introduce cr/l by Cf);*:) = e'’‘‘"^CR{x) and 
cx(.r) = then the Hamiltonian becomes 

Hip) = Vf X P^lo-zCp, (2) 

p 

where Cp = {cRp,CLpY is the Fourier transformation of 
cr/l(x), and cr^ is the spin operator. This Hamiltonian de¬ 
scribes a one-dimensional massless Dirac metal (more pre¬ 
cisely, a Weyl metal), with cr^ playing the role of chirality. 

In this paper we focus on the possibility of dynamical 
fermion mass generation and symmetry breaking ||4^l4^ . For 
simplicity let us take the interaction to be short-ranged, 
namely. Hi - -g(c]^CR)(c^c/,) (with g > 0, which means re¬ 
pulsive interaction among electrons with opposite chirality). 
In the mean field theory, we define m — gic^^CR), and obtain 
its mean field value as \m\ - mo = VFAexp(-2^), where 
A is a momentum cutoff. In this mean field calculation, only 
|m| can be obtained, while the phase of m is arbitrary. It is 
convenient to write m{x) - niQ exp(i0(x)), neglecting the fluc¬ 
tuation of |m|. The ground state without Goldstone modes is 
0(x) = 00. where Oq is an arbitrary real constant. 

The expectation of the x-component of spin is given by 

<cr^(x)> = e^'^''<4(x)cR(x)> -t- e'^''<4(x)cL(x)> 

= 2—cos(2x-0(x)) (3) 

8 

where Q = 2kf. Similarly, we have 

<cr„(x)) = -!e“'®''<4(x)c«(x)) + ie‘^'‘{cl(x)cL(x)} 








2 


s/ 




V ^ A ' 

FIG. 1: Helical order at the edge of QSH. 


= -2—sin(2.t:-6l(x)) (4) 

g 

and {(T^ix)) - {c^{x)cr^c{x)) - 0. The helical order is illus¬ 
trated in Fig. O- So far a physically intuitive mean-field anal¬ 
ysis is presented. A more rigorous approach is the bosoniza- 
tion. The helical spin ordering is favored in the presence 
of sufficiently strong interaction and certain perturbation en¬ 
abling the umklapp processes (see the online Supplemental 
Material). 

Helical order at the surface of weak TI. A natural recipe 
to avoid strong quantum fluctuations in one-dimensional (ID) 
is to couple many ID systems together to form a 2D system. 
This picture brings us to the study of this section. 

Weak TIs are characterized by the so-called weak topolog¬ 
ical indices llsil Is^l . The simplest models of weak topologi¬ 
cal insulators consist of layered QSH [see FiglSJ. Therefore, 
the surface states can be obtained by coupling the QSH edge 
states. Suppose that the surface of the weak TI coincides with 
the xz-plane. To simplify the problem, we include hopping 
terms among only adjacent layers. The low energy Hamilto¬ 
nian for the surface reads 

i/(k) = Vpikx - kF)cj,^qk - Vpjkx + kf)c],^C|k 

k k 

~2f|| cos kj(C||^ci-ii + ^Ik^ik)> (5) 
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FIG. 2: (a) Weak TI as layered QSH. (b) Fermi surface nesting at 
the surface (in the xz plane) of weak TI. 
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where k s (kx, kf), and fy is the interlayer hopping. 

A notable feature of this Hamiltonian is as follows. The 
energy of electrons with spin up and spin down is £’|(k) = 
vpikx - kf) - 2f|| cos kj and £’j,(k) = -Vf{kx + kf) - 2t\\ cos k^, 
respectively, therefore, we have 

Ftlk -1 Q) = Vfikx + kf) + 2f|| cos k^ = -£'i(k) (6) 

where Q s (2kf, n). Therefore, there is perfect Fermi sur¬ 
face nesting at wavevector Q (This salient feature is absent in 
the surface states of the strong topological insulator lfssll l. An 
infinitesimal interaction Hj - -g(cj^CR)(c^cz,) can generate a 
gap, analogous to the case of the QSH edge. The order pa¬ 
rameter of this symmetry breaking is m(x) = g{c]^(x)cR(x)) s 
rriQ exp(!0(x)). The spin densities are given by 

<cr;,(x)) = e"''2'’‘<4(x)cfi(x)) -H e'‘^ *<4(x)cz,(x)) 

= 2— cos(Q ■ X - 0(x)) (7) 

8 

where Q = (2kf, n), and 

<cr^(x)) = -!e"'‘^’‘<4(x)cR(x)) + ie'^'^(cl(x)cL(x)) 

- -2— sin(Q ■ X - 0(x)) (8) 

8 


FIG. 3: Helical order at the surface of weak TIs. This spin order is 
helical in the x direction and staggered in the z direction. 


The spin texture at the weak TI surface is illustrated in Fig[3] 
(We have taken 0(x) to be a constant). The spin order is helical 
in the x direction, while staggered in the z direction. Finally, 
we remark that in real materials the perfect nesting is replaced 
by approximate nesting. To be more conclusive, we have stud¬ 
ied a realistic lattice model in the random phase approxima¬ 
tion (RPA), and found that the approximate nesting favors he¬ 
lical spin order even if the interaction is quite weak(see Sup¬ 
plemental Material), therefore, we expect it to occur in real 
materials of the weak topological insulator. 

Helical spin order in Weyl semimetals. The dynamical mass 
generation induced by interaction and the resultant charge 
density wave state has been studied before I ^ ^ 54-571. 
With an external magnetic field, the Fermi surface instabil¬ 
ity becomes infinitesimal, i.e. an infinitesimal interaction can 
open up a gap at the Fermi surface lfisl 581. The charge density 
wave pattern is, however, too crude to identify the unique fea¬ 
ture of Weyl semimetals, namely, the spin-momentum locking 
described by the Weyl equation. Here we show the existence 
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of helical spin orders, which provides a finer signature of the 
Weyl-type spectrum. 

For simplicity let us consider a single pair of Weyl points 
located at Ki and K 2 respectively. We define the shorthand 
notation Q = Ki - K 2 . The low-energy Weyl Hamiltonian 
reads hAp) = iii which i = 1,2, p,- s 

ki - Ksj, and is a 3 x 3 matrix. For each i = x,y,z we 
can define a vector e* = ep; thus hs becomes more 

compact; 

hs(p) = t'F ^ Pii^l ■ o')> (9) 

i-x,y,z 

where p, is small compared to |Q|. Hereafter we shall focus 
on the cases with e* • for which compact analytical 

treatment is possible. 

Now we add a magnetic field B = Bz (We can always 
rotate the coordinate system such that B points in the z 
direction) 13 • The energy eigenvalues for nonzero Landau 
levels are 

Eiip) = +Vf ^p2 + 2eBl (/ = 1,2,3,...), (10) 

which are all gapped. To obtain the zeroth Landau level, first 
we can solve the 2D problem in the xy-plane by letting p^ = 0. 
In the Landau gauge the 2D Hamiltonian is 

hioiqx, qy) = vfHpx + eBy)e* ■ cr + pj,e* ■ o-)]. (11) 

We can find that the zeroth Landau level wavefunction 
^Px(x,y) = -^= exp[- 5 ^(eBy -H p,t)^] exp(ip^x)|ej x ej), 
where we have introduced the notation |n) to denote the two- 
component spinor satisfying (n|cr,|n) = n, for an vector n. 
Adding the p^ term is now straightforward because |ej x ej) is 
an eigenvector of p^e^ ■ <t (This is the simplification of taking 
■ ej = 6ij). The single-particle wavefunciton is 

X exp(ip;f* + ipzZ)\XseD, (12) 

in which we have introduced the chirality 

= (e:xe;)-e,^ = +l, (13) 

and the single-particle energy is 

Kx.Pz^XsVfPz, (14) 

According to the above wavefunction structure, the fermion 
operators can be expanded as 

c(x) = Yj e'^''’‘[L^e,")c.(x) h- ..., (15) 

S 

where Cj(x) are low-energy fermion operators analogous to 
cr/l(x) in the QSH section, and “.. .” denotes high energy 
modes. Suppose that xi - ~X 2 - T then the index 


identification 1(2) «-> R{L) is valid, and the analysis of 
dynamical symmetry breaking in the QSH section applies, 
namely, a four-fermion interaction induces a chiral conden¬ 
sation (c|(x)c^(x)) = = Y exp(/0(x)). The x-component 

of spin density becomes 

(o-xOi)) = <C+(x)cr^c(x)> 

= e^'^'’‘'^^<X 2 el\o-x\xiel} + H.c. (16) 

s 

If we write \x\^\) - le^) = [cos(0i/2), sin(0i/2)e''^']^ 
and = I - Cj) = [sin((^2/2), - cos((^2/2)e'^"]^, 

then we have {X 2 ez\crx\xi^l} - sin(0i/2)sin((^2/2)e'^‘ - 
cos(0i/2) cos((^2/2)e“'^L Similarly, we have 

(o-yCx)) = <C^(x)cryC(x)> 

= e^'^'^’^^^(X2el\o-y\xiel} + H.C., (17) 

s 

with (xie^lo-ylxi^l) = -/[cos(0i/2) cos(02/2)e~''"2 + 

sin(0i/2) sin(02/2)e''^'], and 

<o-z(x)) = {c\x)o-^c(x)} 

= e^‘^'''^^^{X2el\o-^\Xiel) + H.C., (18) 

8 

with = cos(0i/2) sin(02/2) + 

sin(0i/2)008(02/2)^^^*^^“*^^^ Finally, the charge density 
is 

<o-o(x)) = <c^(x)c(x)) 

= + H.C., (19) 

s 

with Ct'zCjUrie') = cos(0i/2)sin(02/2) - 

sin(0i/2) cos(02/2)e'^'“‘“'^^^ 

Studying some special cases helps us to understand these 
results. For instance, we consider 

/ll(p) = Vf{PxO-x- PyO-z+ PzO-yA 

hlip) = VFiPxO-x + PyO-z + PzO-y). ( 20 ) 

It is readily seen that e* = = (0,1,0), and the previous 

general results tell us that 

{a-Ax)} = — sin(Q ■ X - fl(x)), 

8 

(o-yix)} = 0, 

{(Tzix)} = — cos(Q ■ X - 0(x)). (21) 

8 

This case is shown in Fig.®). If we take a different Weyl 
Hamiltonian 

hl(p) = VfCPxO-x + PyO-y + PzCTz), 

hiCp) = Vf(PxO-x - PyO-y + PzO-z), (22) 
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FIG. 4: Helical spin orders in Weyl semimetals (with dynamical 
mass generation) for (a) Hamiltonian in Eg. 120b and (b) Hamiltonian 
in EQ.l l22t . Here Q is taken to be in the z direction. 


the a simple calculation yields 


cos(Q-x-6'(x)), 

8 

<o-v(x)> = sin(Q-x-6'(x)), 

8 

<o'z(x)) = 0. 


(23) 


This case is shown in Fig.(|4j)). The charge density (cro(x)) = 0 
for these two cases. Ea. (l20l i and Ea. (l22l i qualitatively re¬ 
semble the Weyl semimetal materials and models. For in¬ 
stance, the simple model ifisll /t(k) = 2g sin A:^crj.-i-[2fv(cos ky- 
cos ko) + m(2 - cos - cos k^)]<Ty + 2t^ sin k^cr^ has a pair of 
Weyl nodes at (0, ±ko, 0), with hs=i ^2 = VxCr^Px ± VyCTyPy + 
^z'^zPz its low-energy approximation, which is the same as 
Ea. (l22l l except for the possible velocity anisotropy. A more 
quantitative study of these materials shall be presented else¬ 
where. 

By changing the direction of the magnetic field, |ep is 
changed accordingly, and the helical spin texture changes in a 
prescribed way. We also remark that if the Pauli matrices in 
Eq.(0 are not associated with spin but some other degrees of 
freedom, the helical order is a straightforward generalization 
of the above results. 

The spin helix predicted here can be observed by the spin- 
resolved scanning tunneling microscope (STM). It is unclear 
whether the magnitude of electron-electron interaction (and 
the sample quality) of the recently discovered Weyl semimet¬ 
als favors the generation of the spin helix, but we are probably 
justified in being optimistic about its possible realization, con¬ 
sidering the ongoing rapid progress in this field. 

Electric manipulation of spin texture. So far we have 
not investigated the effects of fluctuations of 0. Such 
phase fluctuations are termed “axions”l5^^, and have also 
been studied in the context of topological insulators and 


superconductors . In our present study, the axions are 

much more visible because of their simple geometrical mean¬ 
ing: they are the phase angle of spin polarization (rotated from 
Q x). 

Now we shall show that spin textures carry electric charge; 
moreover, the charge density depends on the spin texture in a 
precise manner. For concreteness, let us take the QSH edge as 
an example. In the presence of 0(x, f), the spin polarization is 
pointing to angle Q ■ x - 0(x, t). Let us consider the simplest 
case, 6 - 6q-A cos(^x), where A << 1 is the amplitude of the 
spin modulation on the background of the helical order. For 
such a slow modulation of the phase 9 of the Dirac mass, the 
Goldstone-Wilczek formula llhSi] relates the gradient of 9 to the 
charge density: 

1 A 

p(x) ^ —dx9 ^ —‘I sin(?x). (24) 

2n In 

On the other hand, if the charge density is given as p{x) - 
po cos(qx), then we have 

0(x) = —— sin(q'x). (25) 

q 

Therefore, gating the system periodically, such that the charge 
density modulates periodically, can control the spin modula¬ 
tion in a predicable manner. 

We can also consider gating the system to induce a con¬ 
stant charge density po(x) = C. According to the Goldstone- 
Wilczek formula, we have 6 (x) = 27rCx for this case. Now 
the phase Qx — 9(x) = (Q — 2nC)x, namely, the wave vector of 
helical spin order becomes Q - 2nC. This is consistent with 
the relation Q = 2kf: a constant charge density amounts to 
shifting kf in the underlying Fermi surface “before” dynami¬ 
cal mass generation. 

Finally, we remark that taking kf = 0 brings us back to the 
result of Ref. Iholl ; namely, a magnetic domain wall between 
+x and -X magnetization generates a fractional charge ±e/2. 
For a general kf + 0, the +x magnetization is replaced by a 
spin helix with a wavevector Q - 2kf, i.e. both sides of the 
domain wall are spin helices, with a phase difference n. 

Conclusions. A most prominent feature of topological 
Dirac and Weyl semimetals is the spin-momentum locking, 
which is a dramatic consequence of spin-orbit coupling. We 
have shown that this spin-momentum locking can be “frozen” 
as helical spin ordering in the presence of dynamical insta¬ 
bility (“helical solids” from helical liquids). These real-space 
(not the reciprocal-space) helical spin textures should be visi¬ 
ble in the spin-resolved STM. Apart from its intrinsic interest, 
it has potential applications due to the electric tunability of 
helical spin texture. 
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